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Abstract 

Between early 1997 and late 2002, the HEGRA collaboration operated a stereoscopic system of 4 (later 5) imaging 
atmospheric Cherenkov telescopes. In this paper we present the calibration schemes which were developed for the 
system, and report on the performance of the detector over the years. In general, the telescope system was very 
well understood, regarding both the absolute calibration and the slight changes in performance over the years. 
The system had an energy threshold of 500 GeV for observations at zenith and under optimum detector conditions. 
With the corresponding calibration schemes, a systematic accuracy of 15 percent on the absolute energy scale has 
been achieved. The continuous sensitivity monitoring provided a relative accuracy of a few percent, and showed 
that the threshold did not exceed 600 GeV throughout the entire operation time. The readout electronics and the 
imaging quality of the dishes were well monitored and stable. The absolute pointing had an accuracy of at least 
25"; this number was guaranteed throughout the whole lifetime of the experiment. 

Key words: Imaging air Cherenkov technique, Very high energy gamma ray astronomy 
PACS: 95.55.Ka, 95.55.Vj, 96.40.Pq 
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1. Introduction 

The HEGRA 1 stereoscopic Imaging Atmo- 
spheric Cherenkov Telescope (IACT) system was 
located in the Canary Islands, at 2200 m above 
sea level on the Roque de los Muchachos on La 
Palma (17°52'34" West, 28°45'34" North). It 
consisted of 5 identical telescopes (CT2 - CT6), 
which operated in coincidence for the stereoscopic 
detection of air showers in the atmosphere. The 
HEGRA collaboration also ran a single telescope 
(CT 1 |1I2| ) in stand-alone mode, which is not 
covered in this paper. 

The atmospheric Cherenkov technique uses the 
Cherenkov light emitted by the charged shower 
particles, in order to reconstruct arrival direction 
and energy of the primary TeV photons that ini- 
tiated the showers. At the same time, a large 
background consisting of showers induced by the 
isotropically incoming charged cosmic ray parti- 
cles is detected. Since the TeV photons point 
back to their sources, information on the sites 
which emit TeV radiation is directly obtained. 

HEGRA has published a large number of de- 
tailed analyses of characteristics of Galactic and 
extragalactic sources, some of them containing 
data sets which were collected over several years. 
Analysis techniques and the performance of the 
IACT system as well as Monte Carlo simula- 
tions have been described in detail in |3l4l5lfi| . 
Now, after the disassembly of the system, we felt 
compelled to provide a detailed account of the 
long-term performance and stability of the system 
throughout several years of operation. Data for 
this topic have been collected mostly via the stan- 
dard calibration procedure, some information was 
gathered using special calibration setups. Earlier 
results, based on data which were recorded until 
the end of 2000, were in parts already presented 
in|Zj. 

The calibration of the telescopes comprises the 
following major issues: 

Continuous camera electronics calibra- 
tion: The relative gain of the camera pixels and 
the pixel timing was calibrated by so called laser 
runs, in which the whole camera was illuminated 

1 HEGRA stands for "High Energy Gamma Ray Astron- 
omy" 



uniformly by laser flashes. These runs were per- 
formed every night. 

Pointing calibration: The telescopes' point- 
ing was calibrated and monitored by dedicated 
calibration runs, so called point runs, which used 
stars as reference objects. Those calibration runs 
were typically performed every few months. 

Imaging properties: The quality of the op- 
tical point spread function (psf) of the dishes was 
monitored using the same point runs. The influ- 
ence of the psf on 7-ray induced shower images 
could best be obtained from regular observations 
of strong TeV 7-ray sources. 

Relative sensitivity monitoring: Continu- 
ous monitoring of changes in the light sensitivity 
of the complete detectors, including mirror reflec- 
tivity and quantum efficiency of the PMs, could 
be obtained by the measurement of the event trig- 
ger rate which originates from the steady flux of 
cosmic rays. Additionally, the absolute gain of 
the camera electronics alone was monitored by 
the same laser runs as mentioned above. 

Absolute energy threshold: The absolute 
light sensitivity of the telescopes, and hence the 
energy threshold of the system, was addressed 
by so called muon runs which were performed 
roughly once or twice per year, and also by special 
laser setups using a calibrated reference photodi- 
ode. 

These topics are covered in more detail in this 
paper, after an overview on the telescopes' hard- 
ware and an introduction to the data processing 
and general analysis techniques in the next two 
sections. 

2. The telescope system hardware 

Figure shows one system telescope and a 
front view of its camera. Each telescope has an 
altitude- azimuth mount. The reflector has a di- 
ameter of 3.40 m, with a focal length of 4.92 m. 
The mirror is segmented, each dish holds 30 
spherical glass mirrors with diameters of 60 cm; 
the total mirror area per telescope is 8.5m 2 . 
The mirrors are aluminized and quartz-coated. 
The mirror segments are arranged in a so called 
Davies-Cotton design [5], resulting in a design 
point spread function having a spot size of better 
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Figure 1. Up: One of the system telescopes. The 
readout electronics is located in a container which 
can be seen in the background of the telescope. 
Down: Front view of the camera; a funnel plate 
in front of the hexagonally arranged pixel matrix 
collects the incoming light onto the photocath- 
odes. 



than 3' (width of a two-dimensional Gaussian dis- 
tribution) on the optical axis and not worse than 
4' at the edge of the field of view. 

The telescopes are driven by stepper motors 
with a step size of l'/3. The absolute position of 




current monitoring branch T 



Figure 2. Schematic overview of the processing of 
the electric signals from a camera pixel. 



the telescope axes are measured using 14-bit shaft 
encoders (1.'3 per bit). The time used for tracking 
is provided by the local telescope CPUs which are 
synchronized to a GPS clock. The accuracy with 
which the CPU clocks are synchronized is one sec- 
ond. The online tracking accuracy is better than 
0?1 in most cases, being limited by the mechanical 
distortions of the telescope's structures. After ap- 
plication of offline corrections which account for 
these effects, the pointing accuracy is better than 
25" 0. 

The camera in the focal plane consists of 271 
pixels, which are arranged in a hexagonal close 
packed layout. Each pixel contains a photomulti- 
plier (PM) (EMI 9083 KFLA) and a preamplifier 
located in the base. A metallized funnel plate in 
the focal plane concentrates the incoming photons 
onto the photocathodes; the funnels are shaped 
similar to Winston cones. The distance between 
the centers of two adjacent pixels is ~ 21 mm, cor- 
responding to 0?2445 at the sky. The total field of 
view covers a hexagon with side length of 20 cm, 
equivalent in area to a circle of 4?3 diameter. Ca- 
bles provide high voltage (HV) supply from the 
electronics to the camera and signal transfer back. 
For each single pixel, the HV can be adjusted sep- 
arately within a range of 285 V below a camera- 
wide maximum HV value; the pixel values across 
all cameras range from 600 to 1000 V. 

Both AC and DC content of the PM signal 
are amplified separately in the camera, and then 
transferred to the readout electronics (see Fig-EJ • 
which is housed by a container a few meters away 
from each telescope. Except for the HV sup- 
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ply and the next neighbour trigger module, all 
electronics are built in VME standard; the VME 
bus is controlled by a 68040 Motorola CPU. The 
fast pulses induced by air shower events are digi- 
tized with 120MHz, 8-bit flash ADCs (FADC); 
the pulses are shaped with a time constant of 
t = 12 ns in order to match the time resolution 
of the FADCs. The DC currents of the pixels 
are measured for monitoring purposes: During 
data taking, pixels with a DC current above 3 /iA 
(basically because of stars traversing the field of 
view) are excluded from the trigger decision to 
avoid an excessive level of random triggers. These 
pixels are also excluded from offline image anal- 
ysis. Besides that, the current measurement is 
useful for the pointing calibration (see §[5J. 

The trigger decision for an event is made in 
two stages E3 : A local trigger is released when 
two adjacent pixels exceed a certain threshold. 
The comparator thresholds were initially set to 
10 mV, later on to 8mV, corresponding to 8 and 
6 photoelectrons (ph.e.), respectively. The "next- 
neighbour" -decision is done hard- wired with a co- 
incidence window of 12 ns; the hardware solu- 
tion provides fast triggers and allows a low pixel 
threshold. The cameras thus trigger with a rate 
of the order of 10 3 — 10 4 Hz. All of these local trig- 
gers are sent to a central station. If at least two 
telescopes have fired within 70 ns (corrected for 
the zenith angle dependent propagation delay), 
an event trigger is sent back to all telescopes, 
with a frequency of typically 15 Hz. Here, the 
continuous signal digitization by the FADCs is 
stopped at all telescopes, including those without 
a local trigger. For each event, the FADCs and 
the information on which pixels have triggered 
are read out. From the FADC buffers, a window 
of 16 bins per channel (corresponding to 128 ns) 
containing the event pulse is extracted. After a 
zero-suppression, the event data are then sent to 
the central processing computer which performs 
online event building and stores the data to disk. 
The event time is determined from a GPS time 
receiver, with an accuracy of 200 ns; for check 
reasons and possibly as fall-back, also the data 
from a local rubidium clock is read out and stored 
|11| . The dead time after an event is 200 /zs (non- 
extendable, meaning that further events trigger- 



ing during that readout time do not prolong the 
dead time interval). 

The first four system telescopes (CT 3 - CT 6) 
were installed during 1995 and 1996 The sys- 
tem was running in the final hardware layout be- 
tween March 1997 and September 2002. Major 
changes during that time were the reduction of 
the pixel threshold from 8 photoelectrons (ph.e.) 
to the final value of 6 ph.e. in May 1997, and the 
upgrade of the CT 2 electronics and its inclusion 
into the system as fifth telescope in August 1998; 
however, technical problems with CT2 allowed 
the operation of a complete stable 5-telescope sys- 
tem only since June 1999. 



3. Data processing, analysis techniques 
and performance 

During data taking, various preliminary anal- 
ysis steps were performed, including an online 
event display and a fast search for 7-ray signals 
after each data run of typically 20 or 30 minutes. 
For final processing, the data were brought via 
tape from La Palma to a home institute at the 
end of an observation period (3 weeks). In the 
following, the processing of a normal data run is 
described: 

3.1. Pulse shape analysis (pixel level) 

For each event, the FADC pulses have to be 
transformed into pixel amplitudes. First, the zero 
offsets of the FADCs (pedestals) of each pixel are 
subtracted; the pedestals are determined dynam- 
ically with a frequency of a few seconds, from 
the first two FADC bins which, by design, do 
not contain any shower information. A deconvo- 
lution function is applied to restore the original 
pulse. Then the time jitter between individual 
pixels (few ns) is corrected, based on calibration 
results which are discussed in §0] An event time 
is determined from all active pixels above a cer- 
tain threshold. Finally, pixel amplitudes c V i X are 
computed from two FADC bins around that event 
time. Pixels with overflow values in the FADCs 
are treated with a different procedure. Details 
are described in ^21 an d EH- 
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3.2. Image processing (single telescope 
level) 

Pixel processing: The pixel amplitudes are 
transformed into calibrated values, in units 
of the number of photoelectrons (ph.e.) 
which the PM's photocathode has emit- 
ted. For that purpose, all pixel amplitudes 
are multiplied with their respective rela- 
tive gain factors, which are determined from 
the preceding laser run (camera flatficld- 
ing), and a camera- wide conversion factor 
K e i; the latter is determined on a monthly 
basis. Also, a non-linear response function 
for high amplitudes above ~ 150 ph.e. is ap- 
plied. Pixels which are defective or have a 
DC current above a certain threshold are 
set to zero. Details of this procedure are 
discussed in §0] (camera electronics) and in 
S I7.1I (conversion factor). 

Image parameters: To distinguish the shower 
image in the camera from noise pixels, a 
so called two-level tail cut is applied: all 
pixels above 6 ph.e., and all pixels above 
3 ph.e. which have an adjacent pixel above 
6 ph.e. are accepted. From the pixel ampli- 
tudes, the following Hillas parameters ^3] 
are computed for each telescope: center of 
gravity of the image x, orientation of the 
major image axis $, width w, length I and 
the image amplitude amp (sometimes called 
image size). Also, the amplitudes of the 
two "hottest" pixels are stored, i.e. the pix- 
els with the highest amplitudes which have 
(most probably) caused the local telescope 
trigger. These values can be used e.g. to 
track the nonlinear behaviour of the read- 
out chain at high amplitudes, or to apply a 
software threshold similar to the hardware 
trigger condition (cf. § 17. 4|) . 

An event data record finally consists of the im- 
age parameters of all telescopes, the event time, 
and the current actual orientation of each tele- 
scope; these values are written to a DST ("Data 
Summary Table") file. For most analysis pur- 
poses, only DST data are required, which provide 
a fast and easy access to the data. 



3.3. Shower reconstruction (combining 
data from all telescopes) 

From the image data, the shower geometry, the 
shape of the shower, and the energy of the pri- 
mary particle are reconstructed. Here, we con- 
centrate on the following standard analysis chain 
in conjunction with its relevant calibration issues: 

3.3.1. Direction and core location 

The shower axis in space is obtained from a 
pure stereoscopic reconstruction; different algo- 
rithms are described in . The direction of the 
incoming primary particle, for 7-rays equivalent 
to the 7-ray source location, is directly obtained, 
since at TeV energies flight direction of the pri- 
mary and the shower axis are practically identi- 
cal. 

Before the shower axis is reconstructed by any 
stereoscopic algorithm, the positions of the image 
centers Xtei are corrected for the mispointings of 
the individual telescopes. Details of this proce- 
dure are discussed in §[SJ After these corrections, 
all cameras view exactly the same sky sector. The 
shower direction is computed directly in this sky 
coordinate frame, onto which the shower images 
are projected by the telescope imaging, by deter- 
mining the point in the sky where all projected 
shower axes are crossing. 

The reconstructed shower directions are then 
used to produce a sky map. For point sources, 
usually the squared angular distance 2 of the 
shower direction to the source location in the 
field of view (FoV) is plotted. In this repre- 
sentation, a pure background distribution is ex- 
pected to be flat, while 7-ray sources produce 
an excess at zero. All events below an angular 
cut # m ax are counted as ON source events. The 
background estimate is usually derived from other 
areas A in the FoV, and the normalisation fac- 
tor aoN/OFF = ^4on/^4off reflects the different 
opening angles for ON and OFF source counting 
in the FoV. 

The intersection point of the shower axis with 
the ground (i.e. the projected impact point of the 
shower on ground) is called the shower core. The 
shower core is computed similarly to the shower 
direction, in a projection of the images onto a 
plane at ground level which is parallel to the sys- 
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Figure 3. Distribution of the image parameter 
mean scaled width: The left, light-shaded his- 
togram shows 7-rays, derived from background- 
subtracted on-source observations on Mkn501. 
The distribution can be described by a Gaussian 
with a mean of 1.0 and a width of 0.1; a better fit 
can be obtained by a double-Gaussian (solid line). 
The right, dark-shaded histogram shows back- 
ground events, i.e. mostly charged cosmic rays, 
derived from control positions in the FoV; the 
distribution was normalized to the on-source area 
which was used for the 7-ray events. 



tern pointing. For advanced energy reconstruc- 
tion algorithms, the height of the shower maxi- 
mum is determined in addition |16j . 

The reconstructed core locations are used in 
the subsequent shape and energy calculations. 
Note that because of the stereoscopic view of a 
shower with a telescope system, the geometric re- 
construction does not depend on the results of 
those following analysis steps. 

3.3.2. Primary particle identification by 
the shape parameter mean scaled 
width 

Particle identification - in most cases on sta- 
tistical grounds - is provided by the difference 
of the shower development between 7-rays and 
charged cosmic rays (CR), which induce a per- 
manent isotropic background. Usually, a sim- 
ple quantifier like the so called mean scaled width 
(msw) 0] is defined to provide information how 
7-ray-like an event was. 

The mean scaled width is defined as the average 



mSW = > SWtel 

Ntel tr 

tel 



with 



SWtel 



Wtel 



w eKp (amptei, dtctei, zen) ' 



(1) 



(2) 



where w cxp is the expected value for w in a 7- 
ray event, given the image amplitude amp te i, the 
core distance from the telescope dtctei, and the 
zenith angle zen. The sum runs over all N te i tele- 
scopes which have an amplitude amp te i greater 
than 40ph.e. To enhance a 7-ray signal, events 
with a msw value above a certain threshold (typ- 
ically 1.1 or 1.2) are rejected (see Fig.|HJl, the cor- 
responding cut efficiencies are given in Tabled 

The expectation values w cxp were initially de- 
rived from Monte Carlo simulations [H], which 
perform full shower and detector simulations. 
Monitoring and fine-tuning of the values for w exp 
are discussed in §0 

3.3.3. Energy reconstruction 

An energy estimate, presuming that the event 
was induced by a 7-ray, is similarly computed 
from measured event parameters and expectation 
values derived from Monte Carlos [4l5l6ll6lT7| . 

In the standard approach, at first an energy 
estimate E te i is derived for each individual tele- 
scope. The values are derived from a lookup table 



Etei = f(ampl el ,dtc t ei,zeri), 



(3) 



where dtctei is the core distance, zen the zenith 
angle, and amp* tel is a rescaled value of the mea- 
sured image amplitude amp te i ■ The values for the 
lookup table were derived from Monte Carlo sim- 
ulations. The simulations were adjusted accord- 
ing to the status of the experiment in summer 
1997. Subsequent changes in the light sensitiv- 
ity of the detectors are accounted for by rescaling 
the measured amplitude values amp te i by a fac- 
tor «o pt ; the rescaled values are labeled amp* tel . 
More details on the procedure are given in §0 

The energy estimate E of an event is then cal- 
culated as 



E 



(4) 
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The sum runs over all N te i telescopes which 
have an amplitude ampL greater than 40ph.e.* 
and are not too far away from the core location 
(dtctei < 200m for zenith angles below 50°). 

For the spectral evaluation, only events are al- 
lowed for which at least two telescopes are in- 
cluded in the energy estimate, and where at least 
two telescopes have seen the core under angles 
which differ by more than 20°. With the given 
parameters, a 20% single event energy resolution 
is achieved, almost independent of the primary 
7-ray energy |4I5I17|. 

3.3.4. Spectral evaluation 

To derive the source energy spectrum, the 
events are appropriately weighted and filled into 
a histogram. The weight function contains two 
pieces of information: (a) the so-called effective 
area A e s(E), which reflects the energy dependent 
acceptance function of the instrument, and (b) 
the high level analysis cut efficiency k 7 (-E), which 
is imposed by the mean scaled width cut and the 
angular cut. The differential photon flux per en- 
ergy bin Ei is then derived according to 0] 

r\(b 1 f JV ° N,i 

JVoFF.i "I 

-"on/off Vh( E i) A <&{ E i)\~ X ( ■ ( 5 ) 

For each energy bin, the sums run over the 
events in the on-source and off-source areas in 
the FoV, respectively. 

The effective area A c g takes into account the 
pure trigger acceptance of the instrument A' oS [5], 
calculated as 

/>oo 

A' eS (E) = 2w P- l (E,r)rdr (6) 
Jo 

where P 7 (E,r) is the trigger probability for 7- 
rays with a given energy E and a core distance r 
from the system center. Moreover, also the finite 
energy resolution and low level analysis cuts need 
to be considered. This is done by integrating over 
the energy resolution function p(E, E), which de- 
notes the probability that a true shower energy E 
is reconstructed as E and includes the low level 



cuts (basically the cuts on ampl el and on dtc te u 
see § 13.3. 3[1 . A cS is then derived according to 0] 

A M a-,E) = l^^m^l. (7) 

$o, is an assumed model spectrum; for many pur- 
poses the exact shape of this spectrum is not im- 
portant. <I> Q is characterized here by a spectral 
index a of an assumed power law; however, $ a 
can be any model function. Usually one starts 
with a power law with a spectral index a (e.g. 
2.5), and the final spectrum is derived iteratively 
according to the measurement jS] . 

The effective area A e g needs to be calculated 
separately for different system setups (number of 
telescopes actually running in the system), for dif- 
ferent trigger settings, and also for different trig- 
ger sensitivities due to gain changes of the system. 
The latter is done by a simple rescaling procedure; 
details are discussed in §0 

3.4. Performance summary 

In Tabled we summarize the properties of the 
HEGRA telescope system; the values are given 
for a 4-telescope system setup and observations 
at zenith. To first order, the main difference 
between the final 5- and the 4-telescope system 
comes from an increase in the detection rate of 
about 15%; since the HEGRA system has run for 
a considerable time as a 4-telescope-system, most 
calculations and values refer to this setup. 

With the HEGRA system, observations be- 
tween 0° and 60° zenith angle were performed. 
A detailed description of the performance depen- 
dence on zenith angle is beyond the scope of this 
paper. For performance monitoring, experimen- 
tal data up to different maximum angles are used 
in this paper, the limits were chosen according to 
the degree of the respective zenith angle depen- 
dence. Generally, the CR background event rate 
changes only by a few percent between 0° and 
30°. The 7-ray detection rate drops however by 
~ 20% over this range, 7-ray rate comparisons 
therefore need a finer zenith angle binning (0°- 
20°, 20°-30°). 
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Table 1 

HEGRA IACT system performance for a 4 tele- 
scope setup, observations at zenith, and un- 
der optimum detector conditions. FoV: field 
of view. BG: background. t: single event 
resolution, sigma of a 2-dimcnsional Gaussian 
exp (—(9% + #y)/(2cr 2 )). ft : single event resolu- 
tion in one ground coordinate. ftf: single event 
resolution, AE/E nearly independent of E. The 
cut efficiency for 7-rays is called k 7 , kqr is the 
respective fraction of background events after the 
cut. 



year 


mjd 


periods 


runs 


1997 


50449- 


56-68 


5195-8931 


1998 


50814- 


69-80 


8932-13417 


1999 


51179- 


81-92 


13418-17661 


2000 


51544- 


93-105 


17662-21600 


2001 


51910- 


106-117 


22000-27005 


2002 


52334- 


118-126 


27006-29928 



Table 2 

Throughout the Figures of this paper, the time 
development of detector parameters is given in 
one of the above listed units. Each observing pe- 
riod comprises three weeks around the new moon 
date. 



4. Continuous camera electronics calibra- 
tion 

For data calibration and processing purposes, 
monitor and status information about the hard- 
ware of the cameras is obtained mainly from laser 
calibration runs and from normal data runs. 

First of all, pixel channels which are defective 
or found impossible to calibrate are identified and 
are on a run by run basis artificially set to zero. 
Not more than 2% of all 271 channels per tele- 
scope were excluded in most (95%) of the runs. 
In addition, pixels which had in the last moni- 
toring status information a DC current of more 
than 3 fiA (due to stars in the FoV) are also set 
to zero. The DC current values are provided to- 
gether with the single pixel trigger rates in the 
normal data stream; the information per pixel is 
updated online every 16 seconds. 

We note that the exclusion of single faulty pix- 
els does not impose strong inhomogeneity prob- 
lems in the reconstructed sky maps. Since defec- 
tive pixels generally do not appear at the same 
sky positions for different telescopes, possible ef- 
fects are averaged out. Moreover, source candi- 
dates were generally observed offset from the FoV 
center by 0?5, resulting in further averaging of 
camera inhomogeneities because of field rotation 
(due to the telescope alt-az mount). On the other 
hand, if pixels are set to zero due to stars, these 
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Figure 4. Results of a pixel timing calibration for 
one camera. Upper panel: relative arrival times of 
the signals in the FADC, derived from one laser 
run. Middle panel: change of the arrival times 
due to a typical HV flatficlding adjustment for all 
camera pixels. Lower panel: change of the arrival 
times due to an exchange of two FADC modules; 
since the pixels in the camera are numbered in 
concentric circles, not all 16 channels per module 
are adjacent in the corresponding channel plot. 



averaging effects do not apply; noticeable distor- 
tions in the reconstructed sky maps of up to 20% 
in the case of bright stars can occur. This issue 
will be addressed in a forthcoming paper. 

The exclusion of whole telescopes is deferred 
until the shower reconstruction is performed on 



DST level. Here, single telescopes with techni- 
cal problems are excluded. The most common 
reason for that is the number of defective pixels. 
Telescopes with generally more than 15 defective 
channels are excluded. In this case, usually at 
least one of the electronics modules - each com- 
prising 16 channels - is broken. Since pixels which 
belong to the same module are clustered in the 
camera, the loss of one module leads to an inac- 
ceptable inhomogencity in the respective camera. 

The calibration procedure of the PMs and the 
respective electronics readout which follows is 
mainly based on the results of laser runs which 
are performed at the beginning of each observ- 
ing night and occasionally repeated later on. For 
these calibration runs, nitrogen laser pulses are 
fed at each telescope into a scintillator, which 
in return illuminates the camera homogeneously 
with a spectrum which is roughly similar to the 
Cherenkov light spectrum and time distribution 
of real showers. Each laser run comprises 100 
laser flashes with a duration of a few nanoseconds 
and typical amplitudes of 80 to 100 ph. e. These 
values are far above the pedestal values, but be- 
low the level of 150 ph.e. where non-linear effects 
begin to play a role. 

From the laser runs, the relative gains of the 
pixels are directly determined. These results are 
used first of all to adjust every few months the 
HV of all individual PMs. The aim is to provide 
during data taking a trigger acceptance as homo- 
geneous as possible across the camera surfaces; 
this is needed to achieve good sensitivity for the 
compact shower images. This HV fiatfielding pro- 
cedure leaves a 5% trigger acceptance variation 
across the camera FoVs, due to sensitivity differ- 
ences between the FADCs (which are calibrated) 
and discriminators (which provide the trigger de- 
cision). Again, the field rotation and the averag- 
ing over many cameras, in conjunction with the 
usual software thresholds of 40 ph.e. per image, 
led to the fact that no strong inhomogencities in 
sky maps were induced by the trigger threshold 
differences. 

The gain differences between pixels in the 
FADC readout which remain despite the HV ad- 
justment are corrected during data calibration, 
on the basis of the latest laser run result. The 
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pixel amplitudes are multiplied with their respec- 
tive correction factors, which are close to unity 
with a RMS of about 5%. 

Also from the laser runs, the arrival time dif- 
ferences between pixels (a few ns) can be deter- 
mined, since the light from the laser arrives simul- 
taneously at all PMs (Fig.|3| upper panel). The 
signal propagation time depends on the HV value 
of a pixel, resulting in different PM transit times 
(cf. Fig. 01 middle panel), as well as on the hard- 
ware settings of the respective FADC channel (cf. 
Fig-El lower panel). These time differences are 
used in the pulse shape analysis which calculates 
the pixel amplitudes. The following procedure is 
applied: (a) in a first step, an event time is de- 
rived, using only pixels with a reasonably strong 
signal; (b) in a second step, the respective ex- 
pected time at which the signal is expected in 
each pixel is calculated from the event time. In 
both cases, the relative propagation delay is cor- 
rected. In this way, only the true event time is 
sampled also at channels with very weak pulses, 
alleviating the influence of the night sky noise. 

The pixel timing constants were usually derived 
on a monthly basis. After HV adjustments or ex- 
changes of PMs or electronics modules, new ta- 
bles were produced as well. The timing calibra- 
tion even permitted measurement of time profiles 
of air showers, with a duration of a few nanosec- 
onds in the cameras |13| . 

This calibration scheme has worked robustly 
throughout the entire operational lifespan of the 
system. 

5. Telescope pointing 

During data taking, the tracking algorithm 
which steers the telescopes relies solely on the 
position of the telescopes' axes, which are mea- 
sured by optical shaft encoders. The steering 
implementation is able to keep the telescopes on 
track within ±1 shaft encoder unit (0?022). How- 
ever, inaccuracies of the telescopes' mechanical 
structures as well as non-linearities of the shaft 
encoders lead to a typical online mispointing of 
« 0?05, sometimes ranging up to w 0?3. These 
pointing deviations are only depending on the ori- 
entations of the telescopes' axes, and are fully re- 
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Figure 5. Evolution of two pointing model pa- 
rameters for CT3 and CT4. The vertical dotted 
lines indicate the time spans during which one 
respective set of parameters was valid. Within 
errors, the vertical axis tilt has not changed for 
any telescope (dashed line). Shaft encoder off- 
sets have changed quite a few times due to repair 
work; additionally, for CT 4 a small drift of the 
offset was observed after 1999. 



producible. Hence, they can be calibrated and 
corrected offline. 

Since the corrections are small compared to the 
camera size, online corrections are not required. 
The coincidence rate loss and the system FoV in- 
homogeneities, caused by changing overlaps of the 
camera FoVs at their borders, are also small. 

5.1. Pointing corrections 

The pointing of the system telescopes is cor- 
rected using an analytical model of the tele- 
scopes' mechanical structure. This model allows 
one to compute the mispointing of each telescope 
for any given elevation and azimuth angle. It 
parametrizes the bending of the masts and the 
telescope structure, zero offsets of the optical axis 
and of the shaft encoders, tilts of the vertical and 
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the horizontal axes, and the first harmonic devi- 
ation from linearity of the shaft encoders. For 
each telescope, 11 parameters fully describe the 
pointing behaviour. The same model works for all 
system telescopes, but different model parameters 
apply for each of the five telescopes. The system- 
atic pointing error which remains after correction 
is about 25" in both right ascension and declina- 
tion; this value was derived from the calibration 
procedure itself. 

The procedure to derive the model parameters 
makes use of observing selected bright stars as ref- 
erence sources with the telescopes, by performing 
so-called point runs. In one point run, a star spot 
is scanned across the central pixel of the camera; 
13 horizontal scan lines cover a 0?8 x 0?8 window. 
The DC currents of the single pixels are used to 
determine the actual pointing relative to the star. 
One calibration procedure comprises a sample of 
stars distributed homogeneously across the sky. 
More details are given in [5]. 

The initial pointing calibration was completed 
during the telescopes' commissioning phase. Af- 
terwards, the model parameters needed to be re- 
determined only on special occasions, such as a 
telescope mirror readjustment or a dismounting 
and remounting of the shaft encoders. Figure 
shows the temporal evolution of two exemplary 
parameters, the vertical axis tilt and the zero off- 
set of the azimuth shaft encoder, for two tele- 
scopes over their full lifetime. The vertical axis 
tilts of all telescopes are of the order of 0?05, in 
agreement with the accuracy which was envisaged 
during installation, and have, within errors, not 
changed at any telescope. The shaft encoder zero 
offsets, on the other hand, have changed a few 
times due to repair work. For CT 4, an addi- 
tional small drift of the offset was observed (see 
Fig-El j presumably starting after the last manual 
adjustment of this shaft encoder in 1999. 

For the shower reconstruction, the mispointing 
is computed for each event, and the image cen- 
ter of gravities xtei (see § I3.3.1|I are corrected in 
the camera plane. In principle, the mispointing 
also leads to errors in the image orientations dtei ■ 
However, the actual mispointing values are small 
enough that the respective errors of litei are be- 
low the accuracy which can be achieved, and the 



errors are therefore not corrected. 



5.2. The TeV center of gravity of the Crab 

The temporal stability and accuracy of the 
pointing calibration throughout the entire oper- 
ating lifetime of the telescope system was verified 
for example by the examination of the center of 
gravity of the TeV emission of the Crab Nebula. 
Figure shows the sky region around the Crab 
Nebula in celestial coordinates. For reference, the 
Chandra X-ray image is shown. The apparent po- 
sition of the Crab pulsar, being in the center of 
the X-ray image, was placed in Fig.Elat the year 
2000-coordinates of the pulsar, as measured in 
the radio band. The center of the TeV emission 
is determined by a fit of a two-dimensional Gaus- 
sian distribution to all events surviving a 7-ray 
cut (msw < 1.1), reconstructed in celestial co- 
ordinates. The typical angular resolution of the 
I ACT system is 0.1°, but with sufficiently high 
event statistics the center of the TeV emission 
can be extracted with much better accuracy. In 
the case of Crab, the statistical error of the recon- 
structed TeV center position is 3.2 arcseconds in 
each coordinate for the total HEGRA data sam- 
ple (white cross in Fig.^J). 

For a check of possible systematic errors in the 
telescope pointing, the HEGRA data were split 
into six observation periods, from the year 1997 
to the year 2002. Figure shows the deviation 
from the expected position in declination (inset 
with horizontal axis) and right ascension (inset 
with vertical axis), as a function of time. From 
the data one can conclude that within the sum of 
statistical errors (« 7" per year, 3'.'2 for the to- 
tal sample) and systematic errors (25"), the cen- 
ters of the X-ray and TeV 7-ray emission are well 
in agreement over the full lifetime of the exper- 
iment. Although the statistics of the full data 
sample would permit, the systematic error does 
not allow to assign the center of TeV emission to 
either the pulsar or the X-ray cloud. However, 
the given pointing accuracy allowed to reliably 
search for extended TeV emission from the Crab 
nebula. For that purpose, only 7-ray events with 
a predicted angular resolution of better than 0?03 
were used [T%] . 
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Figure 6. Positions of the X-ray and TeV 7-ray emission from the Crab Nebula, shown in celestial coor- 
dinates. The grey-leveled image shows the recent Chandra X-ray image (courtesy of NASA/CXC/SAO); 
the apparent position of the pulsar in the center of the image was placed at the year 2000-coordinates of 
the pulsar, as measured in the radio band. For the TeV 7-ray emission, HEGRA data after a shape cut 
msw < 1.1 and below 30° zenith angle were used; data between 30° and 60° yield within their reduced 
statistics consistent results. The white cross denotes the reconstructed center of the TeV emission with 
its statistical error for the entire data sample. Additionally, the deviation from the expected position is 
shown as a function of the observation epoch, both for declination (graph with horizontal axis) and for 
right ascension (graph with vertical axis). 



6. The point spread function of the tele- 
scopes 

The quality of the point spread function of each 
telescope is primarily determined by the align- 
ment accuracy of the individual mirror tiles and 
by the slight deformations of the reflector itself at 
different elevations. The mirror tiles arc spherical 
glass mirrors of 60 cm diameter; they are man- 
ually adjusted and fixed with screws. The ad- 



justment was done as follows: The telescope is 
pointed 14 degrees above horizon to a point light 
source in approx. 900 m distance. For each mir- 
ror, the appropriate distance to the mirror dish 
and the alignment is found by minimizing the 
spot size in the focal point, which lies in the focal 
plane corresponding to the light source. After the 
alignment procedure, the camera is then shifted 
such that the front of the funnel plate is focused 
to 8 km distance. Since the focal distance of the 
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telescope is 4.92 m, the correct focal plane is only 
displaced by 2.5 cm from the focussing plane. 

The HEGRA system telescopes are operated in 
two tracking modes: normal and reverse mode. 
Reverse mode means that the telescope is driven 
from (an arbitrarily defined) normal mode across 
the zenith in the opposite direction, and thus op- 
erating upside down. While changing the azimuth 
direction does not noticeably influence the mount 
structure, changing the zenith angle puts different 
stress on the mount. Since the mirror alignment 
was usually done in normal mode close to the 
horizon, it must be ensured that deformations of 
the telescope mounts do not significantly change 
the psf of the reflector across the whole observing 
range. This was part of the mount specifications, 
and regular observations up to 60° zenith angle 
were performed. 

6.1. Direct monitoring of the point spread 
function 

The psf can be characterized by its on-axis 
width; the additional off-axis aberrations are 
purely of geometrical origin and are not influ- 
enced by changes of the mirror adjustment qual- 
ity or reflector bending. The dependence of the 
psf on the distance from the optical axis was 
measured during the commissioning of the tele- 
scopes, by scanning stars across the entire cam- 
eras. The parametrisation thus derived was com- 
pared via Monte Carlo simulations with full ray- 
tracing simulations; the simulated image param- 
eter distributions were fully compatible |19j . 

In general, the imaging quality of all telescopes 
was stable throughout the years, also with the 
help of mirror readjustments which were per- 
formed roughly once every year. Some problems 
occurred due to icing of the telescopes which has 
lead to distortions of the mirror tiles' alignment 
a few times; this could be cured by subsequent 
mirror readjustments. 

The deterioration of the psf was monitored with 
point runs. The on-axis psf can be measured 
by standard point runs which are used for point- 
ing calibrations (see §|SJ), where a star is scanned 
across the central pixel. Whereas a single point 
run is not sensitive enough to trace single misad- 
justed mirrors, several runs can be combined to 
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Figure 8. The width of the reflector spot size as a 
function of the zenith angle in normal and reverse 
mode, here for CT6 in fall 1998. The spot size is 
characterized by the width of a two-dimensional 
Gaussian with different values for the horizontal 
axis (upper panel) and the vertical axis (lower 
panel). The accuracy of a single measurement 
can be estimated from the scatter; as expected 
from the way of measurement, it is much better 
in horizontal than in vertical direction. The lines 
show fits of empirical representations for the typ- 
ical dish behaviour to the data points. The light 
grey band denotes the range in which experience 
showed that regular operation was guaranteed. 
A distortion of the mirror alignment shifts the 
psf to higher values; if the psf entered the dark 
grey band, a realignment of the mirrors was re- 
commended. 



produce images such as shown in Fig.|7| Here, the 
psf convoluted with the hexagonal surface of the 
central pixel is represented in grey levels; differ- 
ences of the psf at different elevation angles - due 
to the reflector bending - are averaged out. For 
a well-adjusted reflector, on average 20% of all 
photons which should hit the central pixel are - 
due to the finite spot size - imaged into adjacent 
pixels (Fig.0 right panel). In this representation, 
single or several misadjusted mirrors can easily be 
traced, as shown for example in the left panel of 
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Figure 7. Indirect measurement of the point spread function of an individual telescope; both images were 
obtained summing up different sets of point runs. Encoded in grey levels is the average PM current of the 
central pixel as a function of the displacement of the camera center to the respective star position. The 
currents are a measure of the convolution of the psf with the hexagonal surface function. The pixel size 
is indicated by the hexagonal camera structure, the reflector spot size itself is smaller than one pixel as 
indicated by the reconstructed spot size. The numbers in the pixel centers quote the integrated current 
which falls into the respective surface area, normalized to the total distribution; by definition, this is 
equivalent to the percentage of light which hits the respective pixel if a light beam was homogeneously 
distributed across the surface of the central pixel. Hence, for a well-adjusted mirror (right panel), on 
average 80% of the light of a point source (at random position in a pixel) actually hits the respective 
pixel ("Pix 0"), while 20% is spread across the adjacent 6 pixels ("Ring 1"). The presence of misadjusted 
mirrors (as for example 5 mirrors on the left panel) can easily be detected in this representation. 



Fig-El i n this particular case, five mirrors were 
completely misfocused. 

From a single point run, only the width of an 
assumed two-dimensional Gaussian psf can be de- 
termined, for the vertical and the horizontal axis 
of the spot. To derive these numbers, the Gaus- 
sian is convoluted with the hexagonal pixel sur- 
face area and then fitted to the measured DC cur- 
rent in the central pixel. The results of such an 
analysis, derived from all point runs of one com- 
plete pointing calibration campaign, are shown 
in Fig.OH as a function of the zenith angle. Due 
to the way point runs are performed - horizontal 
scan lines of a star spot across the central pixel 
of the camera - the resolution in horizontal direc- 



tion (upper panel of Fig-EJ) is much better than 
in vertical direction (lower panel). 

The best - i.e. smallest - psf is obtained close 
to 90° zenith angle in normal mode, as expected 
since mirror adjustments are done near this posi- 
tion. Yet, as already mentioned above, the mirror 
dishes were designed to be stiff enough to pro- 
vide an adequate psf over the full zenith angle 
range. The measurements showed that the tele- 
scopes CT 3 to CT 6 have a psf with a Gaussian 
width of tTp S f = 0?03..0?04 per axis in the rele- 
vant zenith angle range |20| • CT 2 however has 
a weaker reflector structure; it was designed as a 
prototype for the forthcoming system telescopes, 
and was integrated last into the system in summer 



15 



1998, after having run as a standalone telescope 
for several years before. Then, after initial mea- 
surements, it was decided to align the mirrors in 
vertical telescope position to achieve best align- 
ment quality near typical - i.e. small - observa- 
tion zenith angles. Nevertheless, the psf of CT 2 
has always been inferior to the other telescopes, 
and made a special data treatment necessary (see 

6.2. Shape parameter expectation values 

The quality of the psf is an important parame- 
ter for the shape of the images that are recorded 
from air shower events. Detector simulations have 
been applied to simulated showers in order to pre- 
dict expectation values for image shape parame- 
ters for 7-rays (and also for charged cosmic rays) 
[Sj. The simulated values were extensively com- 
pared to real 7-ray events obtained from the Crab 
nebula, and especially from Mkn 501 during its 
great outburst in summer 1997 which provided a 
quasi background-free 7-ray sample. Simulations 
and experimental values were found in general to 
be in excellent agreement |4lfi| . 

Usually, the information about the shape of 
a shower event is combined into the single pa- 
rameter mean scaled width (msw, see §|3J). Fig- 
ure shows the distribution of this parameter 
for 7-rays and background events, obtained from 
Mkn 501 (background-subtracted) and from off- 
source background measurements, respectively. 
The distribution for 7-rays should, by definition, 
be centered at mswo — 1.0. The shape is roughly 
Gaussian, with a width of a = 0.1; a more ac- 
curate description can be achieved by adding a 
second Gaussian with 20% of the amplitude and 
a width of 02 ~ 0.15. The expected distribution 
for 7-ray events is needed (a) to predict the ac- 
ceptance Kj.msw of the cut which is applied to 
reject cosmic ray events (usually msw < 1.1 or 
msw < 1.2), and (b) in combination with the 
background distribution to calculate optimized 
cuts for the signal search. 

Precise investigations showed that the value 
of mswQ for real 7-ray data without any cor- 
rections actually was 1.033. Hence, in order to 
achieve the highest possible agreement between 
the acceptances predicted from simulations and 
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Figure 9. Upper panel: Time evolution of the 
mean values sw^^ei of the scaled width distri- 
butions. The values were derived from 7-ray 
data from Crab, Mkn 501 and Mkn 421 (cf. lower 
panel). The vertical dotted lines indicate two 
periods with strongly misadjusted mirrors: CT3 
and CT 4 in the first, all telescopes (mainly CT 3) 
in the second case. Lower panel: The mean value 
mswo of the mean scaled width distribution, after 
application of correction factors c msw> tei to indi- 
vidual telescopes, essentially remained constant 
throughout the years. All values in this figure 
were computed using Monte Carlo expectation 
values which already include the correction fac- 
tor c msw = 1-033 (see text for details). 



obtained from real data, the expectation values 
for 7-rays were multiplied by a correction factor 
Cmsw — 1-033; the distributions shown in Fig.|3 
already include this correction. Apart from that, 
the shape expectation values have remained con- 
stant throughout the lifetime of the telescope sys- 
tem and required no further treatment, with the 
exception of (a) a few periods where misfocused 
mirrors led to measurable distortions of the psf 
(see § 16. If) , and (b) a special treatment of CT2. 
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Distortions of the psf (as e.g. shown in Fig.0 
left panel) lead to a loss of events that would oth- 
erwise correctly be identified as 7-rays, since 7- 
rays are expected to produce smaller images than 
charged cosmic rays. Figure shows in the up- 
per panel the mean value swo -te i of the scaled 
width distribution for individual telescopes (i.e. 
before averaging over all telescopes, see §03). The 
values were obtained from 7-ray events from the 
Crab nebula, Mkn501, and Mkn421, in intervals 
of yearly observing seasons with the exception 
of periods with known distortions of the psf at 
some or all telescopes. We note that with the ex- 
ception of short periods and of CT2, the scaled 
widths average within ±2% of 1. Two periods 
with the strongest temporary deviations from the 
expected value 1.0 are indicated in the Figure by 
vertical dotted lines. Besides that, CT2 perma- 
nently shows mismatching values for sw^^ei for 
most of its lifetime. 

System data or data from individual telescopes 
which were affected by distorted point spread 
functions were either discarded, or the respective 
expectation values were corrected by additional 
factors c msw , t ei (e.g. 0]). The values for mswo 
which include all corrections are plotted in Fig. [51 
lower panel. As can be seen, the expectation val- 
ues msvjQ for 7-rays after corrections were essen- 
tially stable at 1.0 throughout the years. Hence, 
a constant 7-ray acceptance n 7tTnsw for the usual 
cuts was provided. However, when including tele- 
scopes with a poorer psf, an increased background 
level due to the inferior separation power of the 
imaging has to be accepted. 



7. Relative sensitivity monitoring 

The overall detector sensitivity needs to be 
monitored continuously in order to perform the 
energy calibration of the instrument as well as to 
calculate the actual collection area for 7-rays. In 
this section, we discuss the calibration relative to 
the reference periods 62-69 (second half of 1997); 
§00 addresses the absolute sensitivity calibration. 



7.1. The photon conversion efficiency: op- 
tical and electronic gain 

The sensitivity of a Cherenkov telescope is 
mainly determined by the total mirror area, the 
reflectivity of the mirrors, the quantum efficiency 
and gain of the PMs, and finally by the trigger 
threshold. We express the efficiency by which the 
Cherenkov photons are converted into signal am- 
plitudes with two values: 

• The optical efficiency K opt comprises the conver- 
sion of the number of Cherenkov photons which 
hit the reflector area under a certain solid angle to 
the number of photoclcctrons released by the PM 
cathode of the corresponding pixel 2 . In this sec- 
tion, we only deal with a relative factor n' opt close 
to 1, which describes the efficiency change rela- 
tive to some nominal value. The values for K,' opt 
are only derived per period on average for all tele- 
scopes in the standard calibration procedure; the 
method is based on the trigger rate development 
of the system as explained in the following sub- 
section f§ l7.2|) . For differences between telescopes 
see the results in 5 18. II 

• The PM plus electronic amplification in terms of 
a photoelectron-to-FADC-count conversion factor 
Kei,t e z can be determined for each individual tele- 
scope. The values can be derived from the laser 
runs discussed in §0] T^. The method uses the 
fact that the variations of the measured pixel am- 
plitudes from laser shot to laser shot - after cor- 
rection of the laser intensity fluctuations - are a 
measure for the average number of released pho- 
toelectrons; other contributions to the width of 
the signal spectrum can be calibrated in the lab- 
oratory. These electronic conversion factors can 
be measured absolutely, and were set close to 
one FADC count /ph. e. during the commissioning 
phase of the telescopes, by adjusting electronic 
gains and HV values of the PMs. The values for 
Kei,tei are determined for the system calibration 
on a period by period basis. 

Since all telescopes show nearly the same time 
dependence of re e i,tei due to common aging pro- 
cesses, an average value K e \ can be reasonably cal- 



2 Ideal imaging is assumed. 
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Figure 10. Upper panel: System trigger rate (converted to a 4-telescope rate), and event rate after a 
software acceptance cut. Each dot represents one data run of typically 20min; runs up to 30° zenith 
angle were used. The lines indicate the ±10% range of the respective expectation values. Middle panel: 
The deduced change of the energy threshold relative to periods 62-69. For the right energy scale, 500 GeV 
was assumed to be the nominal threshold. Lower panel: Photoelectron-to-FADC-count conversion factor 
k i, averaged over all telescopes, and relative optical efficiency K* opt of the system. The vertical dashed 
lines denote global HV increases of all cameras. The star indicates the result of direct mirror reflectivity 
measurements. The plot ranges from January 1997 to September 2002. For details see text. 
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signal 


unit 


x conversion factor 


ps value 


photon flux (300-600 ran) 


[ph.] 


X K opt 


fa 0.12ph.e.*/ph 


nom. photoelectron counts 


[ph.e.*] 


X K opt 


ps 1 ph.e. /ph.e.* 


act. photoelectron counts 


[ph.e.] 


X K e l 


ps lFADC/ph.e. 


FADC counts 


[FADC] 


Kopt — K opt ' K opt 


ps 0.12ph.c./ph 


K tot = K cl • K opt 


ps lFADC/ph.e.* 


K tot = K opt ' K ol 


ps 0.12FADC/ph 



Table 3 

The table sketches how the signal is converted, to- 
gether with the corresponding conversion factors 
and their typical values. In the reconstruction, 
one starts from the bottom beginning with the 
FADC measurement, and computes the number 
of photons to determine the shower energy. 



culated; it is determined according to 

«cl = Ar 1 ; : ( Kcl,CT3 + Y] Kcl.tei ] (8) 

N *-' + 1 V Irf / 

The sum runs over all N te i telescopes included in 
the system; the central telescope CT 3 is assigned 
double weight, which reflects (very roughly) that 
it is overabundant in the events because of its 
position at the system center. 

The values for n e \ and n' opt are both shown in 
Fig-El lower panel, for the whole lifetime of the 
telescope system. With time, k i degrades con- 
tinuously, presumably due to aging of the PM 
dynodes. The rise of K e i at the positions of the 
vertical dashed lines reflect global increases of the 
HV in all cameras; these were performed to com- 
pensate for the previous gain losses. For K opt , a 
continuous slow degradation is observed. From 
the variations of the data between adjacent peri- 
ods, we can conclude that the relative sensitivity 
monitoring has an accuracy of a few percent. 



7.2. The system event rate and the relative 
optical efficiency 

As already mentioned, the optical efficiency 
K opt itself was not measured directly. Instead, the 
cosmic-ray induced event rate is a good measure 
for the overall detector threshold and sensitivity; 
the system practically does not trigger on noise 
events or local muons ^Hj- The event rate R of 
the system is shown in Fig. llQI upper panel; for 
further processing, the rate was already normal- 
ized to a 4-telescope setup rate, using the empir- 
ical relation i?3t i : R±te\ '■ Rstel — 0.8 : 1 : 1.15. 
To determine the detector sensitivity from the 
event rate, this number needs to be corrected 
for weather influences, and the actual dead time 
of the camera/trigger electronics must be taken 
into account. After these (small) corrections, 
any change of the trigger event rate R - rela- 
tive to the nominal rate Rq = 15 Hz - can be 
directly used to recalculate the energy threshold 
-Ethr, relative to the nominal energy threshold of 
E thl>0 = 500 GeV: 

R oc F CR (E > Eh*) °c E- h ]- 7 (9) 

W^thr.O = (i?0/i?f 58 - (10) 

We note that early measurements of the sys- 
tem trigger rate with different hardware threshold 
settings resulted in a power-law behaviour with a 
spectral index of a ss —1.35 for the standard sys- 
tem trigger settings (> 2 neighboured triggered 
pixels in a camera, > 2 telescopes with camera 
trigger) |TJ]]. Monte Carlo simulations showed as 
well trigger rate dependencies with spectral in- 
dices above the value for CRs, albeit on average 
by Aa ~ 0.15 steeper than the experimental data 

Using the value of a = —1.35 for the thresh- 
old estimate in Eq.QJJ would lead to slight dis- 
continuities in the derived development of K opt 
(cf. Fig. E| lower panel). The value of a = —1.7 
(or even a bit lower) is preferred by this inves- 
tigation. To achieve consistency, a threshold de- 
pendence £thr oc K~ 1,25 (cf. Eg . 112(1 would be re- 
quired. Given that the energy threshold does not 
correspond to a step function in the trigger re- 
sponse, one can expect slightly different effective 
spectral indices for the different ways of chang- 
ing the threshold: on the one hand by setting 
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different hardware threshold levels, on the other 
hand by changing the overall gain of the detector. 
Whether this is the reason for the slight discrep- 
ancy in the spectral indices is unknown. However, 
the results of the present analysis are only weakly 
influenced by the exact value of the effective CR 
spectral index. 

The nominal threshold £W,o was calculated by 
simulations for the reference detector conditions 
Kd = lFADC/ph.e. and < pt = 1 ph.e./ph.e.* 0. 
The time development of the energy threshold is 
plotted in Fig. 1 101 middle panel; we can conclude 
that the energy threshold of the telescope system 
has not exceeded 600 GeV throughout the entire 
lifetime of the experiment. 

The detector sensitivity change can be ex- 
plained to a large extent with the change of k c i. 
The remaining correction factor which is needed 
to understand the full change of the sensitivity is 
assigned to the optics (mainly aging of the mir- 
ror surfaces); thereby K r opt is determined [201 ■ By 
definition, the time interval comprising periods 
62-69 has reference detector conditions, 



calculated such that n' tot 62 _69 



^opt 



is 



lFADC/ph.e.* 
on average during this period, with n' tot defined 
as 

K tot = K cl ' K opt- (11) 

We use the simple estimate 

W£thr,o = KotlFADC/ph.c.*])- 1 (12) 

to translate the energy threshold change into an 
efficiency change. 

For the detector calibration, n' opt is computed 
on a monthly basis from the measured rate R and 
the measured electronic conversion factors K e i,te2 
by combining equations 00 E3 IHI and El 



R \ 58 FADC 



\jpt 



,Ro 



I 



Ni 



tel. 



ph.e.* 



— I «el,CT3 + K cl,tei 
V tel j 



(13) 



The absolute scale of n' opt is arbitrary. How- 
ever, an extrapolation back to the start of the 
telescope system (roughly period 47) yields by 
chance 1; n' opt can thus be interpreted as the 



combined quality degradation of all optical com- 
ponents (mirrors, plexiglass cover of the camera, 
funnel plate, and quantum efficiency of the PMs) 
from the initial value. 

Other values which influence the trigger rate 
remained unchanged: the hardware threshold of 
the telescope system has remained unaltered after 
its final adjustment to 8mV in May 1997; also 
the total mirror area and the mirror alignment 
remained basically constant. 

7.3. Energy and spectral acceptance cali- 
bration 

The calibration values K e i,te/ and n' opt are deter- 
mined on a monthly basis and used for the energy 
calibration of the system data in two ways: 
• The image amplitudes are derived according to 



amptei 



amp tel 



— ^ ^ (^cl,ie/) Cpix 



opt/ 



amptei 



(14) 



(15) 



where c P i X are the raw pixel amplitudes and the 
sum runs over all values of (K e i,tel) 1 Cpix above 
the tailcut threshold (see § 13.2(1 . The rescaled val- 
ues amp* are used in the event energy calculation 
(see §E33J. 

• As the efficiency change also leads to a trig- 
ger acceptance variation, the values are used to 
rescale the energy dependence of the trigger prob- 
ability P 7 (E, r), simply by rescaling the effective 
area A e g(E) (see 5 13.3.4(1 according to 

A* s (E) = A cS (K' tot ■ E), (16) 

where A e g is the effective area as derived by 



Monte Carlo simulations (see 5 13.3.4(1 . 

This procedure is valid provided that (a) the 
trigger acceptance functionality Ry(E,r) just 
scales linearly within the experimentally given 
variations of n' tot , and that (b) nonlinear effects 
(such as the tail cut threshold which is for tech- 
nical reasons computed on the ph.e. -scale, see 
Eg . 114(1 do not play a significant role for the given 
range of < pt . 

7.4. Checks on the relative acceptance cal- 
ibration 

To check the derived efficiency values, an arti- 
ficial software threshold was applied to all system 
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Figure 1 1 . Event rates in the Crab field of view 
(< 20° zenith angle), binned in time intervals 
with similar detector performance. Open sym- 
bols denote background (BG) rates, normalized 
to a circle with radius 9 = 0?22. Filled sym- 
bols represent the 7-rate from the position of the 
Crab, with an opening angle of 9 = 0?22; these 
rates are background-subtracted. Squares denote 
the raw rates without any software cut. Triangles 
are rates after a software acceptance threshold, as 
discussed in the text, circles additionally include 
a loose cut of msw < 1.2. To obtain comparable 
results over the years, CT2 was excluded, and 
only data with 4 active telescopes were used. 



data in the following way: (a) the 2nd-brightest 
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Figure 12. Ratios of background energy spectra 
obtained from the Crab FoV. Only 4-telescope 
data (excluding CT 2) up to 25° zenith angle were 
used, a software threshold as discussed in the text 
was applied. The data were split into 5 obser- 
vation seasons 97/98 to 01/02. For each event, 
the energy was reconstructed under the 7-ray hy- 
pothesis. A raw energy spectrum (i.e. without 
spectral acceptance correction) was filled for each 
time interval, and normalized to the total event 
number to correct for different observation times. 
Each of the histograms was then binwise divided 
by the similarly computed histogram of the to- 
tal data sample. Aa denotes the results of power 
law fits oc E Aa to the ratio histograms, shown by 
the dashed lines. Statistical errors of the data are 
below the point size. 



pixel in each camera must exceed 12ph.e.*; the 
unit ph.e.* indicates that the pixel amplitudes 
were divided by n[ ot beforehand, according to 
the image amplitude amp^ el ; (b) at least 2 tele- 
scopes must surpass this camera threshold. By 
this procedure, the electronic trigger decision of 
the telescope system was simulated, albeit with 
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Figure 13. Ratios of 7-ray energy spectra ob- 
tained from the Crab nebula. See the caption 
of Fig.Elf OT how the histograms were obtained. 
Here, events from the direction of the Crab nebula 
were used, after applying a shape cut msw < 1.1 
to enhance the 7-ray signal, and using back- 
ground control regions to subtract the remaining 
CR background. 



a sufficiently higher threshold than the hardware 
value of 6ph.e. This software threshold should 
provide a homogeneous spectral acceptance for 
all recorded data. 

With the events selected in this manner, three 
tests were performed: 

• The event rate of all system data was redeter- 
mined. The result is shown in Fig. llOl uDDer panel 
together with the raw event rate of the system. 
Within the expected accuracy, the rate is now 
constant, as expected. 

• More important, the 7-ray acceptance and 
hence the 7-ray rate for a steady source should 
also be constant. Figure ^2 shows the 7-ray 
rate from the Crab nebula and the corresponding 
background rates for comparison, obtained from 



zenith angles below 20°, with and without soft- 
ware threshold. The Crab nebula is assumed to 
be a steady TeV emitter on these time scales. 
Indeed, the detected 7-ray rate after software 
threshold remained within 10% constant through- 
out the years. 

• Finally, the stability of the spectral acceptance 
of the detector for 7-rays was also checked with 
Crab data. Figures El an d ED show background 
and 7-ray energy spectral ratios from Crab obser- 
vations. The spectra were split up into yearly in- 
tervals, and were then divided binwise by the full 
1997-2002 spectra (neglecting the correlation). 
The histograms were normalized to the respec- 
tive energy-integrated background event ratios to 
correct for the different exposure times. Data up 
to 25° zenith angle were used; the spectral ac- 
ceptance for 7-rays changes quickly within this 
zenith angle range, but the compared time in- 
tervals have nearly the same zenith angle event 
distribution (with the slight exception of period 
01-02) which makes the comparison feasible. A 
fit oc E Aa (dashed lines) was applied to all ratio 
histograms to quantify the deviation from the ex- 
pected flat ratio (results Aa and \ 2 as denoted 
in the figures) , and also a x 2 for a constant value 
(Xaq=o) was computed. Within the given Crab 
7-ray statistics, the spectral 7-ray ratios are flat, 
as expected after the software threshold, and the 
7-ray spectra from different years are statistically 
compatible. For the last period (01-02), a 2a- 
deviation from the expectation is observed; the 
chance probability for the expected flat ratio is 
8%. This deviation is probably caused by the 
different zenith angle distributions of the com- 
pared data, and is presumably not a hint towards 
the start of limitations of the scaling procedure 
for each period. Corresponding studies on Crab 
spectra using the scaled effective areas A* s did 
not show similar deviations from the expectation. 

Also the spectral indices for the background 
spectra agree within Aa = 0.05; however, be- 
cause of the high number of background events, 
the spectra are not fully compatible, either re- 
flecting the different zenith angle distributions 
of the data sets, or indicating the limitations of 
the applied rescaling plus software threshold ap- 
proach. 
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Data after application of the discussed software 
acceptance threshold can be used to directly com- 
pare measured energy spectra of different sources, 
obtained at different time periods and with arbi- 
trary further analysis cuts (e.g. |21p. However, 
in most spectral analyses, the effective areas af- 
ter calibration as discussed above were applied to 
derive energy spectra (e.g. 4 2 2I23| ). 

8. Absolute energy threshold 

To measure 7-ray spectra and fluxes, one needs 
to rely on the absolute energy calibration of the 
detector. In effect, besides shower simulations 
and the atmospheric transmission, the total con- 
version factor between the number of Cherenkov 
photons and the pixel amplitude, 

K tot = K cl ' K opti (17) 

is required for the detector simulation. 

From initial measurements of new detector 
components, the expected value for the total 
efficiency of a telescope was estimated to be 
«tot,tei = 0.12FADC/ph., where K e \ jte i was as- 
sumed to be IFADC/ph.e. Besides uncertainties 
in the derivation of this number, the aging be- 
haviour of the components is not known a priori. 
Hence, an in situ calibration of the instrument is 
of course desirable. 

The standard method is based on the compar- 
ison of the experimentally measured trigger rate 
induced by charged CRs with the predicted rate 
derived from detailed Monte Carlo simulations 
|B| . Due to uncertainties in the CR flux, spectrum 
and composition, this method has an uncertainty 
of « 15% (or maybe even 22% |22|) m the energy 
estimate. 

Within the HEGRA experiment, two alterna- 
tive methods to directly obtain n to t were inves- 
tigated, a muon ring analysis and an installation 
of a stabilized laser using a calibrated photodi- 
ode as reference. The muon ring results are dis- 
cussed in S I8.1I § 18 . 21 summarizes the absolute en- 
ergy threshold investigations. 

8.1. Muon runs 

The Cherenkov radiation which is emitted by 
a muon hitting the telescope reflector or passing 



muon event 




Figure 14. A typical muon event, as seen by a 
camera. Dark pixels have an entry above 3 pho- 
toelectrons after pedestal subtraction; typically, 
the pixel amplitudes do not exceed 15 photoelec- 
trons. 



nearby creates a ring-shaped image in the camera 
(see Fig.^J |25I26| . The total number of pho- 
tons which have hit the reflector can be recon- 
structed from geometrical properties of the ring 
image alone, i.e. from the diameter and position 
of the muon ring and the asymmetry of the ring 
illumination. Hence, the measured intensity of 
the muon ring can be directly used to determine 
the efficiency with which the instrument converts 
photons into pixel amplitudes. 

During normal operation of a telescope system 
like the HEGRA IACT system, only very few 
muons are recorded, since local muons do not re- 
lease a system trigger. Hence, special muon runs 
are required to effectively trigger on muons. For 
HEGRA, a trigger condition of at least 5 pix- 
els, of which at least 2 must be adjacent, above 
6 mV was used for camera trigger, and the sys- 
tem trigger was disabled. The analysis procedure 
is described in |27| . In general, the agreement be- 
tween simulated muon events and data is good. 
However, for a HEGRA-sized reflector, even the 
highest pixel amplitudes in typical muon rings are 
close to the trigger threshold. Therefore, the de- 
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Figure 15. Measurements of the absolute optical 
telescope efficiencies. The data points are results 
from muon runs. In the upper panels, the devel- 
opment of the efficiency vs. time is shown for all 
telescopes individually. Here, the signal ampli- 
tudes were corrected for the electronic conversion 
factor K c \.tei- The fits (solid and dotted lines) de- 
scribe the data very well; details are explained in 
the text. In the lower panel, the system average 
of all telescope efficiencies is plotted, where the 
relative gain drop of the optical efficiency («opt) 
was also corrected. As expected, this corrected 
system average is constant in time. 



termined efficiencies are slightly biased and pa- 
rameter dependent, which can be explained by 
the influence of the night sky background and 



threshold effects. The estimated systematic error 
for the such derived absolute efficiencies Ktot,tei is 
of the order of 10%. 

The comparison between different efficiencies, 
derived both for different telescopes and for dif- 
ferent times, is however accurate to the few per- 
cent level. In Fig. El the telescope efficiencies are 
drawn as a function of time (filled circles). Each 
data point was derived from all muon runs which 
were taken within one or two nights. The statisti- 
cal errors of the efficiencies are below the marker 
size. In the upper part of the figure, the values 
are shown for all individual telescopes. The pixel 
amplitudes were in this case corrected for the in- 
fluence of K e i t tei (and hence expressed in units of 
ph.e.). Therefore, the derived efficiencies are a 
direct measure of the absolute optical efficiencies 

K opt, teZ- 

A common decline of the optical efficiency is 
observed for all telescopes. The solid lines have 
a fixed slope of -3.7% per year, and were only 
normalized to the muon results. This value for 
d^opt/dt was derived on average for the system 
from the trigger rate analysis (see Fig.^31 and 
Ea. l21|) . and describes the results for the individ- 
ual telescopes very well. The difference to fits 
where the slope was allowed to vary freely (dot- 
ted lines) is small, showing that the results of the 
rate development and of the muon ring analysis 
are well consistent. 

Deviations of a few percent of individual data 
points from the average temporal decline are ex- 
pected, as drifts of the electronic conversion fac- 
tors K e i t tel of that order within an observing pe- 
riod are possible. These drifts are not corrected 
as only average values of K e i,tel from within a pe- 
riod are applied to the data. 

As one can also see from Fig.^J the differ- 
ence of the optical efficiencies between telescopes 
is about 5%. These differences have not entered 
the standard calibration procedure. Since in most 
cases at least 4 telescopes were included in the 
system, the impact on the energy scale is negli- 
gible, and the influence on the energy resolution 
is presumably only small for standard energy re- 
constructions. 

In the lower panel of Fig.^J the system average 
of the optical efficiencies K* pt is plotted, derived 
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after correction of the pixel amplitudes according 
to 

c pix — ( K c\,tel ' K opt) c pix (18) 

(compare to § 17.4(1 . As expected, the such cor- 
rected system efficiency is constant in time. The 
solid line is a linear fit with zero slope, which is 
in good agreement with the data and with the fit 
where the slope was allowed to vary (dotted line). 
The system average value 

K* pt [ph.e.7photon] = 10.7% (19) 

is by definition equivalent to the total effi- 
ciency «tot of the reference period, where K' tot — 
1 FADC/ph.e.*. The measured value agrees 
within 15% with the expected value of 12%, which 
was used to perform the detector simulations. 

8.2. Conclusions for the absolute energy 
calibration 

Both methods, the muon ring calibration as 
well as the laser illumination using a calibrated 
photodiode as reference, have the potential to re- 
duce the error on the energy scale down to a few 
percent. The laser method was described in detail 
in [2E]; the result obtained for CT4 at an early 
stage of the system agrees within 10% with the 
muon run result. 

For the HEGRA system however, systematic 
effects limit the evaluation of the parameters in 
the muon ring analysis as well as in the laser cal- 
ibration. From all performed investigations, we 
estimate the uncertainty on the energy scale to 
be of the order of 15%. 

9. The long-term stability and sensitivity 
of the system 

The aging of the HEGRA IACT system com- 
ponents can be estimated from a linear fit to the 
time development of the efficiencies k c \ and n' opt 
(see §0 Fig.HTj|). The main contribution to the 
aging came from a decrease of the PM gain in the 
camera amplification channel of 

dn cl /dt = (-8.0 ± 0.3)% per year; (20) 

this is most probably an aging effect of the PM's 
last dynode. The continuous decrease of K op t 

dK opt /d< = (-3.7 ± 0.3)% per year (21) 
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Figure 16. Sensitivity of the HEGRA system for 
weak sources, derived from Crab data as a func- 
tion of the experiment's lifetime. The values de- 
note the expected significances for a ten hours ob- 
servation, in units of standard deviations a, for a 
7-ray source which has 10% of the Crab's source 
strength. Filled symbols are for zenith angles be- 
low 20°, open symbols for zenith angles between 
20° and 30°. In the upper panel, angular and 
shape cuts were applied as in a standard analysis 
approach for weak sources. In the lower panel, an 
additional software threshold as discussed in the 
text was applied. 



is presumably caused by a deterioration of the re- 
flecting layer of the glass mirrors. This was also 
supported by direct measurements of the reflec- 
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tivity of several mirrors, which were performed in 
August 2000 (Tig. 1101 lower panel, star marker at 
HEGRA observation period 100); these measure- 
ments were performed relative to new mirrors, the 
result is hence expressed in the same units as n' opt . 

The gradient of the time evolution of /t op t was 
also confirmed by the muon run analysis, as dis- 
cussed in §0 The results (filled circles in Fig. 1151 
upper panels) are very well described by the solid 
lines which have a fixed gradient of —3.7% per 
year. 

Finally, in Fig.^D we show the sensitivity of 
the HEGRA telescope system (in a 4-telescope- 
setup), derived from Crab data, as a function of 
the detector lifetime. The sensitivity is expressed 
in units of the expected standard deviation, when 
observing a source with 10% of the Crab's source 
strength for 10 hours. The values were calculated 
with the measured 7-ray and background rates 
from Crab, as shown in Fig. II II In the upper 
panel, standard analysis cuts for the search for 
weak sources were applied. To first order, the 
sensitivity has remained constant. However, as 
can be seen in Fig.^l a slight decrease of the 
sensitivity with time is observed; the fits yield 

da/dt = (-4.1 ± 0.6)% per year. (22) 

The same trend is visible after application of 
the software threshold as discussed in 5 17.41 
which provides a nearly constant 7-ray accep- 
tance (Fig.^J lower panel, compare to Fig. llll) . 
The reason for this behaviour is presumably the 
following: In order to obtain a constant energy 
threshold (within a range of 100 GeV), the elec- 
tronic gain of the detector was increased in the 
long term, to compensate for the losses in mir- 
ror reflectivity (cf. Fig. llO|) . This however led to 
a slight deterioration of the background rejection 
capability and thus to an increasing background 
by charged CR events. 

10. Conclusion 

Throughout the years, the HEGRA IACT sys- 
tem performance was well monitored and found 
to be stable. In the following, we summarize the 
most important technical data of the experiment 
which were found or confirmed by these investi- 



gations: 

• nominal energy threshold: 
E thrfi = 500 GeV ± 15%(sys) 

• actual energy threshold: 

E thr = £ thl , + AE, AE < 100 GeV 

• aging of the electronic amplification chain 
(presumably aging of the PM's last dyn- 
ode): 

dK c i/dt = -(8.0 ± 0.3)% per year 

• aging of the optical components (presum- 
ably mirror aging): 

dn opt /dt = —(3.7 ± 0.3)% per year 

• absolute pointing accuracy: 
ARa < 25 arcsec(sys) 
ADec < 25 arcsec(sys) 

• alignment of the 7-ray expectation value 
for the shape parameter mean scaled width 
with the nominal value of 1: 

Amswo < 2%(sys) 

which results in an uncertainty of the tight 
shape cut efficiency: 
A«7.™m<i.i < l%(sys) 

• flux sensitivity (5 <r, observations at zenith): 
</) min = 0.1 Crab x (i/lOhrs)" 1 / 2 

with a time gradient of about: 
da/dt = (-4.1 ± 0.6)% per year. 

We believe that with the HEGRA system, 
an important contribution has been made to 
the effort of establishing ground based imaging 
Chcrenkov telescopes as precision detectors, well 
suited within the broad range of astronomical in- 
struments. HEGRA has pioneered the stereo- 
scopic aspect of the ground based air Cherenkov 
imaging technique, a basic feature which is now 
being realized in most of the next generation in- 
struments. The experience gained with HEGRA 
will be profitable for the future projects, such as 
the H.E.S.S. telescope system which is currently 
under construction in Namibia. 
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